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Abstract. Grover's algorithm for quantum searching of a database is 
generalized to deal with arbitrary initial amplitude distributions. First 
order linear difference equations are found for the time evolution of the 
amplitudes of the r marked and N — r unmarked states. These equa- 
tions are solved exactly. An expression for the optimal measurement time 
T ~ 0(y/N/r) is derived which is shown to depend only on the initial 
average amplitudes of the marked and unmarked states. A bound on the 
probability of measuring a marked state is derived, which depends only 
on the standard deviation of the initial amplitude distributions of the 
marked or unmarked states. 
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1 Introduction 

The power of Quantum Computation (QC) was most dramatically demonstrated 
in the algorithms of Shor, for the polynomial time solution of the factorization 
problem jlj , and of Grover |2[ || , for a search which can find a marked element in 
an unsorted database of size N, in 0(y/N) steps (compared to O(N) steps on a 
classical computer). The importance of Grover's result stems from the fact that 
it proves the enhanced power of quantum computers compared to classical ones 
for a whole class of problems, for which the bound on the efficiency of classical 
algorithms is known. This is unlike the case of Shor's algorithm, since in spite 
of the fact that no efficient classical algorithm for the factorization problem is 
known, there is no proof that such an algorithm does not exist. 

A large number of related results followed Grover's original paper 0. Among 
these, the efficiency of Grover's algorithm was analyzed and compared to the 
theoretical efficiency limit of quantum computers for such benchmark search 
problems as introduced (before Grover's result by Bennett et al. [||. The 
algorithm was recently shown to be optimal, i.e., to satisfy the theoretical limit 
H . Further developments include the use of Grover's algorithm or slightly mod- 
ified versions of it as the essential step in algorithms that solve a variety of other 
problems such as quantum search for the median B and the minimum |?|] in 
a set of N items, as well as the collision problem |l. It was also shown that 
other search problems which classically require log 2 N evaluations (queries) of 



a black-box function, can be reduced to a single query using Grover's approach 



|9[ |10| . Finally, it was shown that a simple closed formula describes the time 
evolution of the amplitudes of the generalized problem, which includes several 
marked states jo)- As this work is directly relevant to ours, we briefly summarize 
some of its pertinent results. 

Let k(t) [l{t)] denote the amplitude of the marked [unmarked] states in the 
database, r the number of marked states, and u> = 2 arcsin(-\/ V / N) . It was 
shown by Boyer et al. jll| that after t steps of the algorithm, the marked states' 
amplitude increases as: k(t) = sm[a;(t + 1/2)]/ y/r, while at the same time that 
of the unmarked states decreases as: l(t) = cos[w(i + 1/2)]/ ^/W — r. Since N is 
large, the optimal time to measure and complete the calculation is thus after 
T = 0(y/N/r) time steps, when k(t) is maximal. This analysis relies on the 
fact that the initial amplitude distribution is uniform. However, in a variety of 
interesting cases it would be desirable to apply Grover's algorithm to a non- 
uniform initial distribution. Generically, this could arise in situations where the 
search is used as a subroutine in a larger quantum computation, and the input to 
the algorithm can thus not be controlled. Another example would be cases where 
the given initial distribution over the elements is intrinsically non-uniform. 

In this paper we generalize Grover's algorithm to the case in which the initial 
amplitudes are either real or complex and follow any arbitrary distribution. The 
time evolution of the amplitudes is solved exactly for general initial conditions, 
and the efficiency of the algorithm is evaluated. It is found that for generic 
initial conditions, the search algorithm still requires 0(y N/r) steps, with only 
a constant factor compared to the case of a uniform initial distribution jj| . 

The paper is organized as follows. In Sec. || we define the modified Grover 
algorithm and derive difference equations for the time evolution of the amplitudes 
in it. We solve these equations exactly in Sec. ||, and analyze the results in Sec. 
U A summary and conclusions are presented in Sec. ^. 



2 The Recursion Equations 
2.1 The Generalized Algorithm 

Our modified algorithm is essentially Grover's original algorithm, but without 
the initialization step. It thus consists of the following stages: 

1. Use any initial distribution of marked and unmarked states, e.g., the final 
state of any other quantum algorithm (do not initialize the system to the 
uniform distribution) . 

2. Repeat the following steps T times [an expression for T is given in Eq. (|24|)]: 

A. Rotate the marked states by a phase of ir radians. 

B. Rotate all states by n radians around the average amplitude of all states. 
This is done by applying the "inversion about average" operator, repre- 
sented by the following unitary matrix: 



3. Measure the resulting state. 



2.2 The Dynamics 

We will now analyze the time evolution of the amplitudes in the modified algo- 
rithm with a total of N states. Let the marked amplitudes at time t be denoted 
by ki(t), i — l,...,r and the unmarked amplitudes by h{i), i = r + 1, . . . , N, 
where the initial distribution at t — is arbitrary. Without loss of generality we 
assume that the number of marked states satisfies 1 < r < N/2. We denote the 
averages and variances of the amplitudes by: 



marked: k(t) = °k<t) " ^ & 

1 N 1 N 

marked: l(t) = — — £ l 3 {t) af(t) = — — ^ \l 3 {t) - l(t)\ 2 (2) 



N -r ^ Jy ' 1 w N-r 

j—r+l j—r+l 

The key observation is that the entire dynamics dictated by Grover's algorithm 
can be described in full by the time-dependence of the averages. (The variances 
are defined above for convenience, as they are used later in a different context - 



see Section |4.3j .) Formally, let: 

N 



3=1 j=r+l 



l[(N-r)l{t)-rk(t)] . (3) 



C(t) is thus the weighted average over the marked and unmarked states, with 
the minus sign accounting for the it phase difference between them during the 
algorithm iterations. The following theorem then shows that all states evolve 
equally: 

Theorem 1. The time evolution of all amplitudes ( of both marked and unmarked 
states) is independent of the state index, and satisfies: 

h{t + l)=C(t) + ki(t) i = l,...,r (4) 

h{t + l) =C(t)-k(t) i = r + l,...,N (5) 

Proof. - This follows directly from the algorithm. Consider any marked state 
ki(t); this state is flipped to k[(t) — —ki(t), so that the marked average becomes 
k'(t) = i Y^j=i — —k(t). The unmarked states, on the other hand, do not 
flip, so that the total average after the flip is: x(t) = -^[rk'{t) + (N — r) l{t)\ = 
C(t)/2. "Inversion about average" is by definition: k[(t) — > 2x(t) — k[(t) and 



k(t) -> 2x{t) - k{t). Therefore in total: ki(t) -> C(i) + fci(i) and i 4 (t) -> C(i) - 
From this it follows by averaging that: 

Corollary 2. TTie average marked and unmarked amplitudes obey first order 
linear coupled difference equations: 



k(t + l) = C(t)+k(t) (6) 
l(t + l) = C(t)-l(t). (7) 

These equations can be solved for k(t) and l(t), and along with the initial distri- 
bution this yield the exact solution for the dynamics of all amplitudes by using 
Eqs. (|) and (|. 

3 Solution of the Recursion Equations 

The recursion formulae can be solved by a standard diagonalization method for 
arbitrary complex initial conditions. Let: 

y(t) = (k(t),l(t)) , 

and define: 

N-2r 2(N-r) 2r 

The recursion equations @) and ((?]) can be written as: 



v(t + 1) = A • v(t) , A = 



a b 
— c a 



Diagonalization of A yields a solution for v(i), as follows. Let S be the diagonal- 
izing matrix: 

A^^S- 1 AS = ^-° + ), A ± = 7 e ± -. 

Then w(i) = S _1 ■ v(t) satisfies: 

w(i + 1) = A D ■ w(<) , 

with solution: 

w(t) = ((A_) 4 u;_(0),(A + )* W+ (0)) 

where w(0) = O_(0), w+(0)). This yields k(t) and Z(t) from v(t) = S • w(t). 
Diagonalizing A one finds: 

7 = a 2 + be = 1 (8) 

V 

cosu; = a = 1 — 2 — , (9) 



which is identical to the frequency found by Boyer et al. The eigenvectors 
of A are the columns of S: 



1 



1 1 



r 1 
2 \/ N-r 2 



Using this: 



"<-«^ = w(0) = s- • v(0) = t -i^ l(0) + m 



so that: 

v(t) = S 



L^p(o) + ir(o)) e 

This yields finally, after some straightforward algebra: 



N-r 

k(t) = k(0) cos cut + 1(0) \ smut (10) 



l(t) = 1(0) cos tut- k(0)^J ^ sinujt. (11) 

Together with Eqs. (Q) and (|J) this provides the complete exact solution to the 
dynamics of the amplitudes in the generalized Grover algorithm, for arbitrary 
initial conditions. 

4 Analysis 

Next we derive several properties of the amplitudes. 
4.1 Phase Difference 

The averaged amplitudes can be expressed concisely as follows (even when fc(0) 
and Z(0) are complex): 

k(t) = asin(uot + 4>) (12) 
T(t) = j3 cos(ujt + (f>) (13) 

where 



tan0= ^/Z^ ; a 2 =k(0) 2 +T(0f N 



1(0) V N-r' w w r 1 

^=l( f + k(0f-L- (14) 

which shows that there is a it/ 2 phase difference between the marked and un- 
marked amplitudes: when the average marked amplitude is maximal, the average 
unmarked amplitude is minimal, and vice versa. 



4.2 Constant Variance 

Subtracting Eq. (||) from Eq. (||), and subtracting Eq. (||) from Eq. ([?]), one finds: 



ki(t+ 1) - k(t + 1) = ki(t) - k(t) (15) 
+ 1) - l(t + 1) = -[/*(«) - r(t)]. (16) 

This means that: 

Aki = ki(t)-k{t) and Ak = {-!)%{€) -l{t)], (17) 

are constants of motion (time- independent). It follows immediately from the 
definition that the variances a\ and of [cf. Eqs. (|l|) and (Q)] too, are both time- 
independent. 

This allows us to simplify the expression for the time dependence of the 
amplitudes: 



ki(t) = k(t) + Ah (18) 
h(t) = l(t) + {-If Ah, (19) 

where Ak{ and AU are evaluated at t = 0. 

4.3 Maximal Probability of Success and Optimal Number of 
Iterations 

The probability that a marked state will be obtained in the measurement at time 
t at the end of the process is P(t) = Y%=1 l^( f )l 2 - ^ bound on this quantity can 
be derived as follows. Since all the operators used are unitary, the amplitudes 
satisfy the normalization condition: 

r N 

Y,Mt)\ 2 + E i^)i 2 = 1 (2°) 

i— 1 i=r-\-l 



at all times. Using {y — y) 2 = y 2 — y (y is a random variable), we find from 
Eq. ®: 

N N 

E l^)| 2 = (^V-rK 2 + | E h(t)\ 2 /(N-r). 

i=r+l i=r+l 

Let: 

fmax = 1 — (N — r)af, (21) 

a time- independent quantity. Note that in the case of a uniform initial distribu- 
tion of amplitudes of = and P max = 1- Now, P(t) = P max — (N — r)\l(t)\ 2 , so 
that: 



P(t) < P max (22) 

is the required bound. Using the exact solution, we can show that the P m ax 
bound is in fact tight. For, from Eq. (|l3|) it follows that l(T) — when: 



luT + ^=(j + 1/2)tt, 3 = 0,1,2,... (23) 

At these times the bound is reached so that times T satisfying Eq. ( p3| ) are 
optimal for measurement. Note that this conclusion holds only if fc(0)/^(0) is 
real. When fc(0)/Z(0) is complex, the bound is generally not reached since l(t) 
may never vanish. Collecting our results: 

Theorem 3. Given arbitrary initial distributions of r marked and N — r un- 
marked states, with known averages k(0) and 1(0) respectively, fc(0)//(0) real, 
the optimal measurement times are after: 



t _ (j + l/2) 7 r-arctan[^ v /^ 



arccos (l — 2-^-) 



j = 0, 1, 2, ... (24) 



steps, when the probability of obtaining a marked state is P m ax as given by 
Eq. 

An important conclusion is that to determine the optimal measurement 
times, all one needs to know are the average initial amplitudes and the number 
of marked states. The more difficult case when these are unavailable will be con- 
sidered in a separate publication The expansion of Eq. ( p4] ) in r/N <C 1 (at 
j = 0) yields: 

lfc(0) 



T= ~2l?bf + 4V^r- -^/N + (r/N), (25) 

confirming that Grover's algorithm converges in 0(\J N/r) steps for arbitrary 
distributions. The advantage of an initial amplitude distribution with a relatively 



lfc(O) 

which may significantly reduce the required number of steps. 



high average of the marked states is manifested in the constant offset > 



5 Summary and Conclusions 

In this work we generalized Grover's quantum search algorithm to apply for 
initial input distributions which are non-uniform. In fact, it was shown that by 
simply omitting the first step of Grover's original algorithm, wherein a uniform 
superposition is created over all elements in the database, a more general al- 
gorithm results which applies to arbitrary initial distributions. To analyze the 
algorithm, we found that the time evolution of the amplitudes of the marked 
and unmarked states can be described by first-order linear difference equations 
with some special properties. The most important of these is that all amplitudes 



essentially evolve uniformly, with the dynamics being determined completely by 
the average amplitudes. This observation allowed us to find an exact solution for 
the time-evolution of the amplitudes. An important conclusion from this solu- 
tion is that generically the generalized algorithm also has a 0(^jN/r) running 
time, thus being more powerful than any classical algorithm designed to solve 
the same task. 

This work was initiated during the 1997 Elsag-Bailey - I.S.I. Foundation 
research meeting on quantum computation. 
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